Charge 2e Boson Underlies Two - Fluid Model of the Pseudogap in Cuprate 

Superconductors 



On 
O 
O 
(N 

S3 
3 



CO 



I 

C3 
O 

o 



> 

On 
(N 

o 

OS 

o 



- 1— < 



Shiladitya Chakraborty and Philip Phillips 

Department of Physics, University of Illinois 1110 W. Green Street, Urbana, IL 61801, U.S.A. 

(Dated: June 10, 2009) 

Starting from the effective low energy theory of a doped Mott insulatop^^, we show that the 
effective carrier density in the underdoped regime agrees with a two - fluid description. Namely, it has 
distinct temperature independent and thermally activated components. We identify the thermally 
activated component as the bound state of a hole and a charge 2e boson, which occurs naturally 
in the effective theory. The thermally activated unbinding of this state leads to the strange metal 
and subsequent T— linear resistivity. We find that the doping dependence of the binding energy 
is in excellent agreement with the experimentally determined pseudogap energy scale in cuprate 
superconductors. 



The normal state of the high-T c copper oxide super- 
conductors exhibits a variety of anomalous features in 
the underdoped regime which any successful theory of 
these materials must explain. Central to the exotica 
of the underdoped cuprates are the pseudogap^ and 
strange metal phases. These phases are closely linked 
because once the suppression of the density of states at 
the chemical potential, a key experimental signature of 
the pseudogap, ceases at some critical temperature, T*, a 
metallic state ensues. Such behavior is suggestive of a lo- 
calized, or more properly, a 'bound' electronic state that 
is liberated at T*. While the upturn^ of the resistivity 
at low temperatures is consistent with this bound state 
scenario or charge localizatio n 9 ! 10 ' 11 ! 12 a more direct sig- 
nature is the activated temperature dependenc e) 14 ! 15 ' 16 of 
the Hall coefficient. In a Fermi liquid, the inverse of the 
Hall coefficient is a measure of the carrier density which 
of course is independent of temperature. However, in the 
underdoped cuprates, the inverse of the Hall coefficient 
is strongly temperature dependen t 14 ' 15 ' 16 . Gor'kov and 
Teitel'baum 1 ^ observed remarkably that the charge car- 
rier concentration, nuaih extracted from the inverse of 
the Hall coefficient in La 2 _ ;E Sr.,,Cu04 (LSCO) obeys an 
empirical formula, 



nu_ a u(x, T) = n (x) + m{x) exp(-A(x)/T), 



(1) 



appropriate or a two-component or two-fluid system. 
One of the components is independent of temperature, 
no (a;) (a; the doping level) while the other is strongly 
temperature dependent, n\(x) exp(— A(x, T)). The key 
observation here is that the temperature dependence in 
"Hail is carried entirely within A(x,T) which defines 
a characteristic activation energy scale for the system. 
Gor'kov and Teitel'baum's 1 ^ analysis suggests that the 
activation energy is set by the pseudogap energy scale. 
Consequently, the bound component should be liber- 
ated beyond the T* scale for the onset of the pseudo- 
gap. Should "Hail be an accurate representation of the 
effective charge carrier concentration in the cuprates, the 
above observation indicates that the underdoped or pseu- 
dogap phase necessitates a two-fluid description, which 
has been championed 1 ^ recently to explain NMR, inelas- 
tic neutron scattering and thermodynamic measurements 



on these systems. Nonetheless, the microscopic origin of 
the two fluids has not been advanced. That is, there is 
no microscopic prescription for the precise nature of the 
propagating degrees of freedom that underlie the tem- 
perature dependence of nuaii- For example, Gor'kov and 
Teitel'baum 1 ^ attributed the unbinding of the localized 
charges above T* to excitations from van Hove singu- 
larities at the bottom of the band up to the chemical 
potential. 

By contrast, our explanation of the the two fluids re- 
lies entirely on the strong correlations of a doped Mott 
insulator, that is, Mottness. Here we show that the ex- 
act low-energy theory of a doped Mott insulator—^ de- 
scribed by the Hubbard model naturally resolves the two- 
component conundrum in the cuprates. The propagating 
degrees of freedom that constitute the two fluids are the 
standard projected electron in the lower Hubbard band 
and a bound composite excitation composed of a charge 
2e boson and a hole. It is the unbinding of the latter 
that gives rise to the strange metal regime. The binding 
energy is found to be in excellent agreement with exper- 
imental values for the pseudogap energy scale. 

We review some of the key features of the our effective 
low energy theory of Mottness, the complete details of 
which have been worked out elsewherei^i 4 -. Our start- 
ing point is the usual Hubbard model 



fli 



Hubb 



-t 



9ij C i. 



(2) 



where i,j label lattice sites, gtj is equal to one if i,j are 
nearest neighbours, Ci a annihilates an electron with spin 
a on lattice site i, t is the nearest- neighbour hopping ma- 
trix element and U the energy cost when two electrons 
doubly occupy the same site. The cuprates live in the 
strongly coupled regime in which the interactions domi- 
nate as t ~ 0.5eV and U = 4eV. A low-energy effective 
action is then obtained by integrating out the physics 
on the t/-scale. Because double occupancy occurs in 
the ground state, integrating out the [/-scale physics is 
not equivalent to integrating out double occupancy. We 
solved this problem by extending the Hilbert space to 
include a new fermionic oscillator which represents the 
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creation or annihilation of double occpancy only when 
a constraint is solved. The new fermionic oscillator en- 
ters the action with a mass of U and hence represents the 
high-energy scale, which must be integrated out to gener- 
ate the low-energy action. The corresponding low-energy 
theory contains new degrees of freedom, namely a charge 
2e boson, denoted by tpi, that are absent in the original 
model and are not made out of the elemental excitations. 
ipi enters the theory initially as the Lagrange multiplier 
to maintain the constraint that in the extended Hilbert 
space the heavy fermionic field represents the creation of 
double occupancy. To leading order in t /U , the effective 
Hamiltonian 



u 



(3) 



contains the t — J model (the first two terms) and new 
interactions with the charge 2e boson, ipi which represent 
mixing with the sectors with varying numbers of doubly 
occupied sites. 



Here b { = = 



with Vj = —V± = 1. The |6,| 2 term contains the spin- 
spin interaction as well as the three-site hopping term. 
A gradient expansion of this term shows that the spin- 
spin term scales as a 4 , a the lattice constant, whereas the 
terms linear in b are proportional to a 2 . Hence, relative 
to the terms linear in b, the \b\ 2 term is irrelevant. Our 
key contention which has been worked out extensively 
for the cuprates 1 "^ is that as far as the charge degrees 
of freedom are concerned, it is the interactions with the 
tpi sector that determine the propagating degrees of free- 
dom, not the dynamics arising from the spin-spin term. 
In particular, we show that it is the tp terms that give 
rise to a gap in the single-particle electron spectrum. As 
this gap is on the order of t, any contribution from the 
spin-spin term will be subdominant. 

It is important to realise that once the heavy field is 
integrated out, the Hilbert space reverts back that of the 
Hubbard model. Further, as tpi has no bare dynamics and 
tp has no Fock space of its own, its only contribution will 
be to create bound states within the Hilbert space of the 
Hubbard model. This can be seen initially by considering 
how the electron operator transforms 1 * 2 ^^ 



u 



(4) 



upon the integration of the high energy scale. The first 
two terms represent the standard electron operator in 
the lower Hubbard band dressed with spin fluctuations. 
However, the last term represents the correction due to 
dynamical spectral weight transfer—, that is the mix- 
ing with doubly occupied sites. All such processes are 



mediated by ipi which represents a collective charge 2e 
mode arising from the dynamics of double occupancy. 
The quantity p\a- a represents a bound complex with a 
non-propagating local degree of freedom. It is from this 
term that the bound state dynamics emerges. Roughly, 
the two-fluid emerge from the fact that an electron at 
low energies is a linear superposition of an essentially free 
part, the first two terms in Eq. (jJJ and the last term in 
Eq. (3J from which the bound -stated dynamics (that 
is, pseudogap physics) emerges. As a result of tp+, the 
conserved charge is no longer just the total number of 
electrons but 



(5) 



In order to obtain the electron Green function for the 
effective Hamiltonian, we treat <pt as spatially indepen- 
dent, owing to a lack of gradient terms in that field in 
the Hamiltonian. To keep track of the dependence on 
ipi it is helpful to introduce the rescaling tp^ — > sip%. The 
electron Green function is then written as a path integral 
over the ip fields as 

G(k,w) = - J[Dp*][Dp]FT(J[dc*][dc t ] Cl (t)c:(0) 



x exp 



L(c, tp)dt 



(6) 



where the effective Lagrangian L is expressed in a diag- 
onalized form 



klklk, 



(7) 



where the jka are the Boguliubov quasiparticles and are 
given by 



7^ T = cos0fc4 T +sin8 k c^ kl 



in 



sin 6»fcc£ T + cos 9 k c- kl 



(8) 



(9) 



where cos 2 9 k = |(1 + j*-), a k — 2(cosfcj; + cosk v ) , 

E Q = (-2/x + £)<p*(p,E k = -g t ta k - ^,X k = v^ + A 2 , 
the gap is proportional to s, = sip*(l — jj-a k ), and 
hence vanishes when tp is absent and gt = j^pj- S = 1 — n. 
The gt term originates from the correlated hopping term, 
(1 - n lS )c\ cr c jl y{l - rijff). The The 7 fc(T 's play the role of 
the fundamental low energy degrees of freedom in a doped 
Mott insulator. That is, they are the natural propagating 
charge degrees of freedom. Note they depend in a com- 
plicated way on the the tpi field and consequently are 
heavily mixed with the doubly occupied sector. Starting 
from Eq. l[Tj). we integrate over the fermions in Eq. ^ 
to obtain, 
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G(k,u>) = - I [D<p*][D<p}G(k,u,tp)e X p~ 



(10) 



where 



G(k,iv,(p) = 



sin 2 8 k [ 
oj + Xk [ 



cos" 



Xk[p] 



(11) 



is the exact Green function corresponding to the La- 
grangian, Eq. ([7]), which has a two-branch structure, cor- 
responding to the bare electrons and the coupled holon- 
doublon state respectively. The role of the ip field, which 
determines the weight of the second branch, is vital to 
our understanding of the properties of Mott systems, as 
was demonstrated previously-^. It is trivial to see that 
in the limit of vanishing s (no tp field), the 7/ C(T 's reduce 
to the bare electron operators c& and the first term in 
Eq. tfTTj) vanishes. The two-fluid nature of the response 
stems from this fact of the theory. Namely, the first term 
contributes only when (p ^ and the second when (p = 0. 
These contributions correspond to the dynamical and 
static components of the spectral weight, respectively. 

We obtained the Green function G(k, uj) by a numer- 
ical integration of Eq. (fTU|) over the ip field. The Hall 
coefficient Rh was computed from the spectral function 
A(k, oj) using the Kubo formula 19 



t 2 

y XX' 



(12) 



where 



27T 2 \e\ 3 aB f S/H 1 ^,5e k , 



3r 



x !^(k,.) 3 



(13) 



and 



Tre 2 f df(u) s 1 ^£>e k 



KE&M^f (14) 



duo 'N^ y dh 

k 
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FIG. 1: riHaii plotted as a function of inverse temperature 
for four different values of hole doping x: 1) solid circles, 
x — 0.05, 2) diamonds, x — 0.10, 3) triangles, x = 0.15, 
and 4) squares, x = 0.2. The inset shows the temperature 
independent part of the carrier density as a function of x. 
Note it exceeds the nominal doping level indicated by the 
straight line. 
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with <7 XX and cr xy the diagonal and off-diagonal compo- 
nents of the conductivity tensor respectively, f(u>) is the 
Fermi distribution function, and B is the normal compo- 
nent of the external magnetic field. The effective charge 
carrier density riHaXJ is then obtained using the relation 
Rh = -l/("-Haiie). 

FigQ] shows a set of plots of rtHaii as a function of 
the inverse temperature, each corresponding to a differ- 
ent value of hole-doping, x, in the underdoped regime 
(x ranging from 0.05 to 0.20). The plots fit remarkably 
well to an exponentially decaying form. In other words, 
the computed charge carrier density within the charge 
2e boson theory of a doped Mott insulator agrees well 
with the form given in Eq. ([T]) proposed by Gor'kov 



FIG. 2: A(x) (solid circles) obtained from fitting the plots 
in Fig.(fl} to Eq.( [1} plotted as a function of hole doping 
x. The experimental values are also shown for LSCO: solid 
triangle a 15 ' 16 ' 21 and squares^ The excellent agreement indi- 
cates that the bound component contributing to the charge 
density does in fact give rise to the pseudogap. 



and Teitel'baumi^. The inset shows the temperature- 
independent part of the charge density as a function of 
x. This quantity exceeds the nominal doping level. This 
deviation is expected as the Hall coefficient is expected 
to change sign around x = 0.3^ in hole-doped samples. 
The 'binding energy', A(x), was extracted for each 
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FIG. 3: T*(x)) (solid circles) obtained from Eq.([T3J plotted 
as a function of hole doping x. The experimental data were 
gleaned from the following: open circles are from Ref.^T*, 
open triangles (T m ) from Ref.—, and closed triangles {T m ) 
from Ref.—. 

doping and plotted in Fig. (J2J) using Eq. |T]). Shown here 
also are the values for the experimentally determined 
pseudogap energy for LSCQ 15 i 16 ' 21 . The magnitude of 
A(x) falls with increasing hole doping as is seen experi- 
mentally and hence is consistent with its interpretation, 
even quantitatively, as a measure of the pseudogap tem- 
perature T*. A more accurate estimate of T* , 

T*{x) ^ -A(x)/\n(x), (15) 

may be obtained from A(a;), by equating the num- 
ber of doped carriers x with that of the activated 
ones nx(x)exp{— A(x, T)) as proposed by Gor'kov and 
Teitel'baumi-.. Fig.© shows a plot of T* as a function 



of x. This is in reasonable agreement with the experi- 
mentally obtained estimates of T*&22&. 

Ultimately it is not surprising that the pseudogap^- 3 - 
appears within the charge 2e boson theory. As mentioned 
previously, the charge 2e boson is a local collective non- 
propagating mode that is restricted to mediate electronic 
states within the Hilbert space of the Hubbard model. 
The only option is that the boson binds to a hole to form 
a new charge e state. As an electron at low-energies (Eq. 
(HI) is a linear superposition of the standard state in the 
lower Hubbard band and the bound state mediated by 
the charge 2e boson, a two-fluid charge model is a nat- 
ural consequence. This further supports the idea 3 that 
the pseudogap temperature T* represents the boundary 
between bound and unbound charge 2e bosons where the 
binding energy to excite a boson vanishes, and T- linear 
resistivity obtains 3,4 . The mechanism for T— linear re- 
sistivity is simple within this model. Once the binding 
energy of the boson vanishes, bosons arc free to scatter 
off the electrons. The resistivity of electrons scattering 
off of bosons is well-known to scale linearly with temper- 
ature above the energy to create a boson. Hence, this 
mechanism is robust and should persist to high temper- 
atures. Consequently, the charge 2e boson theory offers 
a resolution of the pseudogap and the transition to the 
strange metal regime of the cuprates. 
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